Differentiating the Method of Conjugate Gradients

نویسندگان

  • Serge Gratton
  • David Titley-Péloquin
  • Philippe L. Toint
  • Jean Tshimanga Ilunga
چکیده

The method of conjugate gradients (CG) is widely used for the iterative solution of large sparse systems of equations Ax = b, where A ∈ R is symmetric positive definite. Let xk denote the k–th iterate of CG. This is a nonlinear differentiable function of b. In this paper we obtain expressions for Jk, the Jacobian matrix of xk with respect to b. We use these expressions to obtain bounds on ∥Jk∥2, the spectral norm condition number of xk, and discuss algorithms to compute or estimate Jkv and J T k v for a given vector v.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Conjugate gradient neural network in prediction of clay behavior and parameters sensitivities

The use of artificial neural networks has increased in many areas of engineering. In particular, this method has been applied to many geotechnical engineering problems and demonstrated some degree of success. A review of the literature reveals that it has been used successfully in modeling soil behavior, site characterization, earth retaining structures, settlement of structures, slope stabilit...

متن کامل

Lanczos and Linear Systems Lanczos and Linear Systems

abstract Lanczos's major contributions to the numerical solution of linear equations are contained in two papers: \An Iteration Method for the Solution of the Eigenvalue Problem of Linear Diierential and Integral Operators" and \Solutions of Linear Equations by Minimized Iterations ," the second of which contains the method of conjugate gradients. In this note we retrace Lanczos's journey from ...

متن کامل

Natural Gradients in Practice: Non-Conjugate Variational Inference in Gaussian Process Models

The natural gradient method has been used effectively in conjugate Gaussian process models, but the non-conjugate case has been largely unexplored. We examine how natural gradients can be used in non-conjugate stochastic settings, together with hyperparameter learning. We conclude that the natural gradient can significantly improve performance in terms of wall-clock time. For illconditioned pos...

متن کامل

Parallel Preconditioned Conjugate-Gradients Methods on Transputer Networks

We show how highly efficient parallel implementations of basic linear algebra routines may be used as building blocks to implement efficient higher level algorithms. We discuss the solution of systems of linear equations using a preconditioned Conjugate-Gradients iterative method on a network of transputers. Results are presented for the solution of both dense and sparse systems; the latter bei...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:
  • SIAM J. Matrix Analysis Applications

دوره 35  شماره 

صفحات  -

تاریخ انتشار 2014